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Abstract 

We calculate analytically low frequency quasi-normal modes of gravitational perturba- 
tions of AdS Schwarzschild black holes in d dimensions. We arrive at analytic expressions 
which are in agreement with their counterparts from linearized hydrodynamics in S d ~ 2 x R, 
in accordance with the AdS/CFT correspondence. Our results are also in good agreement 
with results of numerical calculations. 
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1 Introduction 



Quasi-normal modes (QNMs) determine the late-time evolution of black hole perturba- 
tions. They have been extensively studied in asymptotically Anti-de Sitter (AdS) space- 
times [1] in hopes of shedding some light on the Anti-de Sitter - conformal field theory 
(AdS/CFT) correspondence. Analytic expressions for asymptotic frequencies were derived 
in [2, 3] (adapting the monodromy argument proposed in [4] and extended to first order 
in [5]) for arbitrary dimension. First-order corrections to these analytic expressions were 
obtained in [6] in good agreement with numerical results [7]. 

At the other end of the spectrum, the low frequency modes can be used to probe the 
behavior of the gauge theory on the boundary of AdS in the hydrodynamic limit [8]. This 
is particularly interesting in connection with heavy ion collisions at RHIC (see, e.g., [9] and 
references therein) and the LHC. It appears that the quark-gluon plasma (QGP) that forms 
is strongly coupled raising the possibility that it may possess a gravity dual. In [10], the 
QGP was analyzed in terms of a "conformal soliton flow" by mapping the boundary of AdSs, 
S* 3 x R, to the four-dimensional flat Minkowski space by a conformal transformation. The 
QNMs were calculated numerically and led to an elliptic flow coefficient and thermalization 
time that compared well with experimental results. 

Here we calculate the low frequency QNMs of AdS Schwarzschild black holes in d di- 
mensions analytically following the procedure of ref. [8] . We find agreement with numerical 
results in four [7] and five [10] dimensions. Our analytic expressions also agree with their 
counterparts obtained from hydrodynamics on S d ~ 2 x R [10, 11] in accordance with the 
AdS/CFT correspondence. We discuss gravitational perturbations using the Master Equa- 
tion of Ishibashi and Kodama [12], including vector (section 2), scalar (section 3) and tensor 
(section 4) modes. Our conclusions are summarized in section 5. 



2 Vector perturbations 



We wish to study QNMs of AdS Schwarzschild black holes whose metric in d dimensions is 
given by 

ds 2 = -f(r)dt> + ^ + r^ d _ 2 , f(r) = * + 1 - . (1) 

To simplify the notation, we shall set the AdS radius R = 1. The radius of the horizon is 
related to the mass parameter by 

2„ = r*" 1 (l + ^) (2) 

The radial wave equation for gravitational perturbations in the black-hole background (1) 
can be cast into a Schrodinger-like form, 

-^2+V[r{r^ = uH, (3) 

in terms of the tortoise coordinate defined by 

±1 = J_. (4 ) 

dr f(r) { > 
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It is also useful to introduce the distance 



[°° dr 

r* = r*(oo) - r*(0) = J (5) 

A short calculation yields [3] 



(d - l)r H V rf - 1 

The potential V is determined by the type of perturbation and may be deduced from the 
Master Equation of Ishibashi and Kodama [12]. For vector perturbations, we obtain [3] 

V(r) - Vv(r) s m {^LzH + (j - mz mr) _ } (7) 

Evidently, the potential vanishes at the horizon (F(rf/) = 0, since /(rff) = 0). This is the 
case for all types of perturbation. 

The asymptotic form of QNMs was found in [2, 3] analytically. Subsequently, it was 
shown in [6] how the approximation may be improved upon by a perturbative expansion. 
One obtains 

u n = w (p) + w w + . . . 

d-3\ 1, „ 
n+— j.+ -ln2 

, cos4 spfesy (^/r H r^ n^hWH^jd - 3) 

9 2 1 73 — T\ \°) 

27Tr H [2{d-2){l + l/r 2 H )}— (d- 1 ) 

. l{l + d - 3) d 2 - 8d + 13 
^ V= d-2 + 2(2d-15) (9) 



,(0)- 



where 



Thus, the first-order correction is 0(n~^i). It should also be noted that the first-order 
correction is suppressed by a factor of l/rfj. So for large black holes, the zeroth-order 
contribution provides a good approximation to all modes, not just the high overtones. 
Moreover, the zeroth-order term is independent of both the angular momentum quantum 
number I and the type of perturbation. This leads to a mild dependence on I and the type 
of perturbation of all QNMs for large black holes (ru > 1)- Finally, since f* ~ l/rn, the 
frequencies are proportional to the radius of the horizon. 
In particular, for d = 5, the QNMs are given by 

'~ = ( 2n + 1 + i— ) 0(l/r 2 H ) (10) 



rjj \ it 
Numerically, we find 

^«f 3 + i—Vl- = 3.221 -2.779i , ^ = " n+1 - ^ * 2(1 - i) (11) 

These results are independent of the type of perturbation and the quantum number i. They 
agree well with numerical results [10]. 
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Even though the above expressions include high as well as low frequencies, they do not 
include the lowest overtones. The latter may be obtained using the method of [8] and they 
correspond to the hydrodynamic behavior of the gauge theory fluid. 

To obtain analytic expressions for the lowest overtones, it is convenient to introduce the 
coordinate 

..(2)" (12, 

The wave equation (3) becomes 

-(d-3) 2 u^f{u) (u^fiuW)' + V v (u)^ = Co 2 ^ , Cj = — (13) 
V J r H 

where prime denotes differentiation with respect to u and we have defined 

/(.) = ^ = i-.*(.-i£!) (14) 

vvw = * - /<-){* + ™»-/w - «>-m-v( 1 + i/ u } m 

where 

L 2 = ^ + d 2 ~V (16) 

r H 

Let us first consider the large black hole limit r# — > oo keeping u> and L fixed (small). 
Factoring out the behavior at the horizon (u = 1) 

= (1 -u) _i a=iF(u) (17) 

the wave equation simplifies to 

AF" + B^F' + C^ L F = (18) 

where 

„ 2d-8 d-1 

A = -(d-3) 2 u~{l-u~) 

2d-8 d-1 
d-1 d-5 „ iuj U d ~ 3 (1 — U d ~ 3 ) 

= -(d-3)[d-4- (2d-5)u^]u^ -2(d-3) 2 



d-l 1 - it 



f2 (d-2)[d-4-3(d-2)u^l 
C- f = L 2 + ^ ^ — i i rf-a 



" 2 ^ + (d-3) 2 «-"'(!-«--») 



d-1 d-5 2d-8 d-1 

iCj [d - 4 - (2d - 5)u d -z]u d -* 2 iCj m«(1-m«) 

For small u, L, eq. (18) may be solved perturbatively by writing it as 

(H + Hi)F = (19) 
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where 



H F = AF" + B F' + C 0fi F 

HiF = (B G> -B Q )F' + (C Cltl -C Q ,o)F (20) 
Expanding the wavefunction perturbatively, 

F = F + F 1 + ... (21) 



at zeroth order we have 
whose acceptable solution is 



H F = (22) 



d-2 



F = n 2 < d - 3 ) (23) 



d-2 



It is regular at both the horizon (u = 1) and the boundary (u = 0, or ^ ~ r 2 — > as 
r — ► 00). The Wronskian (up to an arbitrary multiplicative constant) is 

W = 1 (24) 

Ud-3 (1 - u<*-3) 



Another linearly independent solution is 



F = F J ^ (25) 

It is unacceptable because it diverges at both the horizon (Fq ~ ln(l — u) for u ~ 1) and 

the boundary (Fq ~ n 2 ( d -3) for u w 0, or ^ ~ r^ - — > 00 as r — > 00). It may be expressed 
in terms of hyper geometric functions but its explicit form is not needed for our purposes 
(first-order perturbation theory). 
At first order we have 

H F l = -HiF (26) 



whose solution may be written as 



Fl - F °J F$j ^4W- (27) 



The limits of the inner integral may be adjusted at will because this amounts to adding an 
arbitrary amount of the unacceptable solution (25). To ensure regularity at the horizon, 
we should choose one of the limits at u = 1 (the integrand is regular at the horizon, by 
design). Then the behavior of the wavefunction (27) at the boundary (u = 0) is given by 

where we omitted regular terms. The coefficient of the singularity ought to vanish, 

- 1 FoT^Fh = 

AW K ' 
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which imposes a constraint on the parameters (dispersion relation) of the form 



a L 2 - iaiuj - a 2 uj 2 = (30) 
After some algebra, we arrive at explicit expressions for the coefficients, 

a = ^— ^ , ai = d - 3 (31) 

The coefficient a 2 may also be found explicitly for each dimension d, but it cannot be 
written as a function of d in closed form. However, it does not contribute to the dispersion 
relation at lowest order. E.g., for d = 4, 5, we obtain, respectively 

a 2 = ^_Il n3 !-Ihi2 (32) 
108 3 ' 6 2 v ' 

The quadratic in Co equation (30) has two solutions, 

1? d — 3 L? 
uJo^-i- — - , wi«-! Vi- — - (33) 

d — 1 3i2 d — 1 

where we omitted terms of higher order in I? . In terms of the frequency lo and the quantum 
number £, they may be written respectively as 

£(£ + d-3) wi d-3 £(^ + d-3) 
(d - l)r H r H a 2 (d - l)r 2 H 

The smaller of the two, uq, is inversely proportional to the radius of the horizon and is 
not included in the spectrum (8) obtained earlier by different means [2,3,6]. The other 
solution, uj\ is a crude estimate of the first overtone in the spectrum (8). Numerically, for 
d = 5, we obtain using (32) 

OL„ -4.109, (35) 

to be compared with the numerical value (11). 

It should be noted that the crude estimate (35) already exhibits two important features: 
u\ is proportional to r# and the dependence on £ is of order as expected from (8). 

The approximation may be improved by including higher-order terms in the perturbative 
expansion. Inclusion of higher orders also increases the degree of the polynomial in the 
dispersion relation (30) whose roots then yield approximate values of more QNMs. Thus, 
this method reproduces the spectrum (8) albeit not in an efficient way. 

The above discussion may be extended to black holes of finite size in a straightforward 
manner by treating finite-size effects as a perturbation (assuming 1/th is small). Thus at 
first order, we need to replace Hi (eq. (20)) by 

Hi = Hi + \h h (36) 

where 

H H F = A H F" + B H F' + C H F (37) 
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The coefficients may be easily deduced by collecting 0(l/r 2 H ) terms in the exact wave 
equation given by (13), (14) and (15). We obtain 



Ch 



-2(d-3) 2 u 2 (l-u) 



-(d-3)u 
d-2 



(d-3)(2-3u) - (d-l) 



1 — U d-1 

U d-S 



d-1 
1 - UdS 



d-A - {2d - 5)u - (d - 1) 



d-1 

1 - u d --i 



U d-1 
U d-3 



(38) 



Interestingly, the zeroth order wavefunction Fq (eq. (23)) is an eigenfunction of Hh, 

H H F = -(d-2)F (39) 
therefore, the first-order finite-size effect is a simple shift of the angular momentum (eq. (16)) 



IS 



L 



d-2 



■ H 



Consequently, the QNMs of lowest frequency (34) are modified to 

.£(£ + d-3) -(d-2) 



CJ 



For d = 4, 5, we have respectively, 



LU 



(d - l)r H 



l)(£ + 2) 



3r H 



+ o(iy H ) 



+ l) 2 -4 
Ar H 



(40) 



(41) 



(42) 



in agreement with numerical results ([7] and [10], respectively). 

Next we discuss the role the lowest frequency mode plays in the AdS/CFT correspon- 
dence. The dual to the AdS Schwarzschild black hole is a gauge theory fluid on the boundary 
of AdS (S d ~ 2 x R). To find the dual of vector perturbations, one ought to consider the fluid 
dynamics ansatz 

u . = JCe~ inT Yi (43) 

where Ui is the (small) velocity of a point in the fluid and Vj is a vector harmonic on S d ~ 2 . 
Demanding that this ansatz satisfy the standard equations of linearized hydrodynamics, 
one arrives at a constraint on the frequency of the perturbation Q which yields [10, 11] 



n = -i 



; £(e + d-3) -(d-2) ( , 
(d-l)r H 



+ 0(l/r 2 H ) 



(44) 



in perfect agreement with its dual counterpart (41). 



3 Scalar perturbations 



Scalar perturbations are also governed by a Schrodinger-like wave equation (3) with the 
potential given by [12] 



V s (r) = 



4 r 2 



£(£ + d-3)-(d-2) + 



(d-l)(d-2)fi 



1 -2 
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d(d - l) 2 (d - 2)V 6(d - l)(d - 2) 2 (d - + d - 3) - (d - 2)]/i 



r.2d-8 



rd— 5 



+(d - 4)(d - 6)[^ + d - 3) - (d - 2)] V + 



2 2{d- l) 2 (d-2)V 



4, ,3 



„3a!-9 



4(d - l)(d - 2){2d 2 - lid + 18)[l(l + d - 3) - (d - 2)}fi 2 

r 2o!-6 

(d - l) 2 (d - 2) 2 (d - 4)(d - 6)/x 2 6(d - 2)(d - 6)[*(£ + d - 3) - (d - 2)] 2 ^ 



„2d-6 



r d— 3 



6(d - l)(d - 2) 2 (d - 4)[^ + d - 3) - (d - 2)]// 



r d— 3 



+4[^(£ + d - 3) - (d - 2)] 3 + d(d - 2)[^ + d - 3) - (d - 2)f 



(45) 



QNMs are asymptotically the same as the QNMs (8) of vector perturbations [2,3]. They 
differ at first order in the perturbative expansion. For scalar modes we obtain [6] 



(i)= - 



'2*3=2) cog 4 7T 



A 5 



d-3 



2ir 2 r 2 H 



— r(— -)r 4 ( ; 



d-3 



where 



^4s 



[2(d - 2)(1 + l/r^)]3=2 * "d - 2'~ v 2(d - 2) * 
(d 2 - 7d + U)[£(£ + d - 3) - (d - 2)] 2d 3 - 24d 2 + 94d - 116 



(d- l)(d-2) 2 



+ 



4(2d-5)(d-2) 



(46) 



(47) 



Again, we see the ^-dependence entering at 0(1/^), as for vector perturbations. Also, the 
above set of frequencies does not exhaust the spectrum as it leaves out the lowest frequency 
mode. To find it, we work as before. Changing variables to (12), we may write the wave 
equation as in (13) with Vv replaced by 



V S (u) 



/(«) 



rh + 1 + -o- ] u 



H 



1 -2 



1 \ 2d-8 ( I \ d-5 

d(d - 2) 1 + -=- J U — - 6(d - 2)(d - 4)m 1 + -5- J ti« 



'H 



' H 



+(d-4)(d-6)m 2 u"^ + (d-2) 2 1 + ^- u 



H 



+2(2d 2 - lld + 18)m ^1 + 
(d-4)(d-6) + 



'H 



U 



u 2 -3(d-2)(d-6)m 2 1 + 



H 



H 



U 



6(d-2)(d-4)m(l + ^) 



m 



u + 2(d-l){d-2)m i + d{d-2)^ r \ (48) 



H 



where 



rh = 2 



£(£ + d-3) - (d-2) 2(i + d-2)(i- 1) 



(d-l)(d-2)r 2 H (d-l)(d-2)r 2 H 
8 



(49) 



In the large black hole limit r# — > oo with m fixed, the potential simplifies to 



d-l ( 

" ln\ 1 — U d ~ 3 2d-S d-5 

r>) = — -,{d{d-2)u^=z - 6{d - 2)(d - 4)mud=s 



4(m + u) 2 
+(d - 4)(d - 6)rh 2 u~^ + (d - 2) 2 u 3 

+2(2d 2 - lid + 18)mu 2 - 3(d - 2)(d - 6)m 2 u + 2(d - l)(d - 2)m 3 I (50) 



In addition to the singularities at the end points (u = 0, 1), the scalar wave equation has 
an additional singularity due to the double pole of the scalar potential at u = —rh. It is 
desirable to factor out the behavior not only at the horizon, but also at the boundary and 
the pole of the scalar potential. We therefore define 



* = {i- u )~^l. F{u) 

K ' m + u K ' 



(51) 



In terms of F(u), the wave equation for scalar perturbations in the large black hole limit 
reads 

AF" + BtjF' + C^F = 



(52) 



where 



A = -(d-3) 2 u^(l-u^) 

2d-8 d-l 

Be = -(d-3)ud-z (l-ud-s 



d - 4 2(d - 3) 



u 



m + u 



2d-8 



d-l . 



d-l d-5 



-(d - 3)[d - 4 - (2d - 5)u^=3}u^=3 - 2(d - 3)' 



iuj u d -3 (1 - ii d -3) 



d-l 



1-u 



Z.U.— Q LI- 

-U d "3 (1 - U d - 



(d-2)(d-4) (d-3)(d-4) 2(d-3) 21 



4u 2 



u(rh + u) (rh + u) 2 

2d-8 d-l 
iCJ U d-3 (1 - U d-3) 



{ % d-l 1 d-5 ^ JJJ-3 M - 
d-4- (2d-5)u— }u~ + 2(d-3)- — 



u 



d-4 d-3 



2u m + u 



-(d - 3) J^ld-*-W-S)u*-hutl _ {d _ 3)2 



2d-8 d-l 
iCj U d - 3 (1 — U d " 3 ) 



d-l 



1-u 



d-l 



(1-u) 2 



1 - u d - :i 



— + (d-3Y 



UJ 2 U d -3 (1 — It d-3) 

(d-l) 2 (1-u) 2 



To calculate the QNMs perturbatively, we define the zeroth-order wave equation as in (22) 
with 

H F = AF" + B Q F' (53) 



The acceptable zeroth-order solution is (arbitrarily normalized) 

F (u) = 1 



(54) 



which is plainly regular at all singular points (u = 0, 1, — rh). On account of (51), it 

d-4 

corresponds to a wavefunction vanishing at the boundary (^ ~ r 2 as r — ► 00, using 
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(12)). The Wronskian is 

W = 

It d-3" (1 - U d -3) 

The unacceptable solution is 

F = j W (56) 
and can be written in terms of hypergeometric functions. It has a singularity at the bound- 

d-5 d-6 

ary, Fq ~ u d -» for it sa (or \£ ~ r 2 — > 00 as r — > 00 for d > 6; for d = 5, the acceptable 
wavefunction behaves as r -1 / 2 whereas the unacceptable wavefunction behaves as r -1 / 2 lnr; 
for d = 4 the roles of Fo and Fo are reversed, however our results are still valid). Fq is also 
singular (logarithmically) at the horizon (u = 1). 

At first order in perturbation theory, working as in the case of vector modes, we arrive 
at the constraint 

(57) 



/ 

J 



/o AW 

where we used (29), (54) and 

H!F = {Be - BoW + C^F = (58) 
The first-order constraint (57) may be written as a dispersion relation 

ao — aiwi — SL2UJ 2 = (59) 

After some algebra, we obtain 

d-l l + (d-2)m d-3 l /1xn ,. u , RfV . 

a = ^ M - M2 , ai= , a 2 = — {1 + O m } 60 

2 (l + m) z (l + m) z m 

For small rh, the quadratic equation (59) yields the solutions 



„± .d — 3 „ d — 1 ^ 

u o ~ m ± y "^J - m ( 61 ) 

The two solutions are related to each other by a)^ = —00^*, which is a general symmetry 
of the spectrum. Neither solution is included in the spectrum obtained from asymptotic 
QNMs. To obtain approximations to those modes, we need to consider higher orders in 
perturbation theory. Notice also that unlike vector modes, these scalar lowest frequency 
modes have finite real part. Using (49), we may express the frequencies in terms of the 
quantum number I as 



± „ _ l(l + d-3)-(d-2) / f(l + d-3)-(d-2) 
J ° ~ Hd 3) (d-l)(d-2)r H V d-2 



Thus the real part is independent of ru whereas the imaginary part is inversely proportional 
to r H (for r H > 1). 

Finite size effects may be added perturbatively. At first order, they amount to a shift 
of the coefficient ao in the dispersion relation (59), 

a — > a + -yag (63) 
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Working as in the vector case, after some tedious but straightforward algebra, we obtain 

aH = ±-{l + 0(m)} (64) 
The modified dispersion relation yields the modes 



- + d — 3 A d — 1 „ 1 

cjq ps — i — - — m± \ — - — m + 1 (65) 



correcting (61). Explicitly, 



± ~ -(a * £(£ + d-3)-(d- 2) / l(l + rf-3) 
"0 "-^" 3) (d-l)(d-2)r g ± V d-2 (66) 

correcting (62). This is in agreement with numerical results [10, 11]. 

Turning to the AdS/CFT correspondence, we ought to purturb the gauge theory fluid 
on the boundary of AdS (S d ~ 2 x R) using the ansatz 

m = JCe- iQr Vi$ , Sp = JC'e- iQT § (67) 

where Ui is the (small) velocity of a point in the fluid, 5p is the pressure perturbation and Sj 
is a scalar harmonic on S d ~ 2 . Demanding that this ansatz satisfy the standard equations of 
linearized hydrodynamics, one arrives at an expression for the frequency of the perturbation 
17 [10, 11] which is in perfect agreement with our analytic result (66). 



4 Tensor perturbations 

Finally we discuss tensor perturbations. In this case the potential is given by [12] 

y Tir) = m { W> + (j - Ml WW + MzgQr)} („) 

It leads to the same asymptotic form of QNMs as in the other two cases (vector and scalar) 
[2,3]. The first-order correction to the asymptotic expression is [6] 

Wr\ [2(d _ 2)(1 + 1/r?;)]3 i i r( < i - 2 )r( 2 ( d - 2 ) ) < 691 

where 

_ t { j + d -3) (d-3) 2 

- 4t " d-2 + 2~(2d^) (70) 

to be compared with the vector (8) and scalar (46) cases. Once again we obtain a mild 
0(l/rfj) dependence on the quantum number I. 

Unlike the other two cases, this constitutes the entire spectrum. To see this, let us 
change variables to (12) and go to the large black hole limit. The wave equation simplifies 
to 

-{d - 3) 2 (u^ - ^ 3 )*" - (d ~ 3)[(d - 4)1*3=1 - (2d - 5)u 2 }^' 

f f2 d{d -2) (d-2) 2 Co 2 \ T 
+ {L 2 + K — > u i-s + —> u — U = (71) 



' 1 -U3=3 
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The zeroth order equation is obtained from (71) by setting L and u> to zero. The resulting 
equation may be solved exactly. Two linearly independent solutions are (* = Fq at zeroth 
order) 

F (u)=u J ^ , F (u) = u~^^\n(l-u^^ (72) 

Neither behaves nicely at both ends (u = 0, 1), therefore both are unacceptable. It is not 
possible to build a perturbation theory to calculate small frequencies. 

This negative result is in agreement with numerical results and is also in accordance with 
the AdS/CFT correspondence [10]. Indeed, there is no ansatz that can be built from tensor 
spherical harmonics Tjj (similar to the vector (43) and scalar (67) cases) satisfying the 
linearized hydrodynamic equations because of the conservation and tracelessness properties 
of Tij [10]. 

5 Conclusion 

We calculated analytically low frequency QNMs of gravitational perturbations of AdS 
Schwarzschild black holes in arbitrary dimension using the Master Equation of Ishibashi 
and Kodama [12]. We noted that low frequency modes are well approximated by asymp- 
totic expressions for large black holes (with radius of horizon r# > 1 in units such that the 
AdS radius R = 1) [2, 3]. These expressions (ujn^ in eq. (8)) are independent of the type of 
perturbation and the angular momentum quantum number I. The dependence on I enters 
at first-order perturbation theory and is 0{\jr\j) [6]. The first-order contribution, u>n \ for 
vector, scalar and tensor perturbations is given by eqs. (8), (46) and (69), respectively. 

Asymptotic expressions do not in general yield the entire spectrum. To find the lowest 
frequency mode, we applied the method in ref. [8]. We also included the effects of a finite 
size black hole. We arrived at explicit analytic expressions in the case of vector (eq. (41)) 
and scalar (eq. (66)) modes. For tensor modes, this method does not yield a new mode 
(the asymptotic series exhausts the spectrum). Our analytic expressions were in agreement 
with numerical results [7, 10, 11]. They also agreed perfectly with the results from linearized 
hydrodynamics of the gauge theory fluid on the boundary of AdS space [10, 11] in accordance 
with the AdS/CFT correspondence. 

It would be interesting to extend these calculations to a less symmetric configuration 
that would better fit the experimental setup of heavy ion collisions at RHIC and the LHC. 
Then a comparison with experimental results would enhance our understanding of gauge 
theory fluid dynamics at strong coupling. 
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